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DEFINITION OF SYMBOLS 


Symbol 

X 


P 


e 

e* 

m 

m* 

c 

a 

M 

A 


CO 


->■ 

CO 


AAA 
X, y, z 


AAA 

T,d,Cp 


M 


Definition 

London penetration depth; the characteristic distance over which currents 
and magnetic fields vary within a superconductor. In the London theory, 
X^ = m* 

P= 1/X; a symbol introduced for mathematical convenience 

coherence length; the characteristic distance over which superconducting 
properties vary within a superconductor 

density (number per unit volume) of superconducting electrons 

electronic charge (note that e < 0) 

charge on a superelectron (Cooper pair); e* = 2e 

mass of the electron 

mass of a superelectron; m* = 2m 

speed of light 

electrical conductivity 

magnetic permeability 

Angstrom unit; a unit of length equal to 10'* ° meter 
angular frequency of rotation 

a vector directed along the spin axis having magnitude cu 

unit vectors in a Cartesian coordinate system (coordinates x, y, z) 

unit vectors in a spherical coordinate system (coordinates r, 0, 0) 

London Moment; the magnetic field outside a rotating superconducting 
sphere is the same as that due to a j^oint dipole located at the center of 
the sphere having a magnetic moment M = (m*c/e*)[ 1 - (3coshj3R/|3R sinhjSR) 

+ (3/f)^R^)]R^Zo— (m*c/e*)R^cS 
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DEFINITION OF SYMBOLS (Concluded) 


Symbol 


Definition 


Pe(x) 


Legendre polynomials; functions orthogonal on the interval [-1,1] 


Ygj^(0,0) spherical harmonics; functions orthonormal on the surface of the unit 
sphere [0 < (^ < 27 t] , [0 < 0 < #] 


P 


electric charge density 


B, H magnetic induction, magnetic field intensity, respectively; related by the 

constituitive relation ^ = /xit 


^ 

e, h, j 


local values of the electric field, magnetic field, and current density, 
respectively. 


N 


torque 


L2 angle between the gyro spin axis and the z-axis of the coordinate system 
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MAGNETIC TORQUE ON A ROTATING SUPERCONDUCTING SPHERE 


I. SUMMARY 


The London theory is reviewed in Section IIIA, and a theorem extremely useful 
for calculating torque on a rotating superconducting sphere is proved in Section lllB. The 
torque on a rotating superconducting sphere due to a uniform applied field is then 
calculated exactly in Section IllC. The rotating superconducting shell is considered in 
Section HID and an estimate of the minimum thickness of the superconducting coating 
for the rotor is obtained. The direct effect of excess charge on the rotor is calculated in 
Section HIE. 

r 

A macroscopic approach is introduced in Section IVA and is shown to give results 
for fields, currents, and torques which are correct to order X/R; the torque in an 
arbitrary magnetic field is then calculated in Section IVB. It is pointed out in Section 
IVC that problems for which the only superconductor present is an infinite 
superconducting plane can be solved by image methods; in Section IVD this fact is used 
to calculate the torque on the rotor due to a current flowing in the readout loop. In 
Section IVE the torque due to a flat spot in the superconducting magnetic shield 
surrounding the rotor is discussed, and the torque is calculated for the limiting case of a 
very large flat spot. The results of the calculations are discussed in Section IVE. 


II. INTRODUCTION 


The scientific goal of the Stanford gyroscope relativity experiment is a test of 
general relativity through measurement of relativistic precessions of Earth-orbiting 
gyroscopes [ 1 ] . The concepts and ideas for implementation of the experiment have been 
described in an earlier document [2]. 

Fundamental design requirements for the experiment are: 

1 . Development of a spin-axis readout with sufficient stability and sensitivity to 
measure the small relativistic precessions. 

2. Reduction of all nonrelativistic torques to less than 2 x 10’*® N m to achieve 
a design goal of less than 1 ms of arc per year residual error. 

Since an asymmetry of some sort is necessary for readout, requirement (2) adds the 
constraint that no significant torque be introduced through the asymmetric feature added 
for this purpose. 



The proposed rotor is a fused quartz ball coated with a thin layer of 
superconducting niobium. The asymmetry for readout is to be the magnetic moment 
generated by a rotating superconductor along its instantaneous spin axis (London 
Moment). Precession will be measured by measuring changes in magnetic flux through a 
superconducting loop encompassing the rotor. 

This report presents results of the author’s calculation of certain torques coupled 
to the rotor through its London Moment. A calculation of the direct effect of excess 
charge on the rotor is also included. 


III. CALCULATIONS BASED ON THE LONDON THEORY 
A. Introduction to the London Theory 

There are two fundamental . lengths associated with superconductivity. The 
characteristic distance over which, currents and magnetic fields vary within a 
superconductor is called the penetration depth, X(T). In the superconducting state, the 
velocities of two electrons are correlated if the distance between them is less than a 
certain range |. This length, ^(T), is the characteristic distance over which superconducting 

properties [for example, the pair potential A(r)] can change and is called the coherence 
length of the superconductor. It should be emphasized that both X and f are temperature 
and material dependent. 

In 1934, F. and H. London introduced a theory of electrodynamics for super- 
conductors in order to explain the Meissner effect [3]. This theory has since been 
shown to follow from the condition of minimum free energy (in the limit of weak fields 
and currents) provided the criterion X » ^ is satisfied [4] . It has further been 
shown [4] that the microscopic (quantum mechanical) theory of superconductivity 
reduces to the Landau-Ginzburg theory for temperatures near the transition temperature 
T^., and the Landau-Ginzburg theory in turn reduces to the London theory for X»^ 

The London theory permits a complete phenomenological description of 
superconductivity which is in accord with the principles of thermodynamics and classical 
electrodynamics. 

The London theory is applicable to the calculation of magnetic torques on the 
gyro: The transition metals generally have large penetration depths (X~ 1000A) and 
small coherence lengths (^ ~ 100 A) at T = 0; niobium is in this class'. Superconducting 
shields around the gyro will ensure that the weak field condition is met. 

1. Actually, for bulk pure niobium, X(0) ~ J(0) ~ 400 A [5-7], However, niobium films 
inevitably have trace impurities which increase X(0) to ~ 1000 A [6,8] , while the shortened 
mean free path decreases ^(0) to ~ 100 A [9]. 
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The London theory assumes Maxwell’s equations^ to be valid inside a 
superconductor: 



— 

47T ^ 

1 dt 

I. 

curl h = 

— j + 

c 

1 ah 

c at 

II. 

curl e = 
— ^ 

•c at 


III. 

div h, = 

0 

> 

IV. 

div ^ = 

47rp 



Here j is the total current density, p is the density of electric charge, and e and h are, 
respectively, the intensities of electric and magnetic fields^. The relationship between 

field intensities and the supercurrent density is postulated to be: 


V. curl (X2 J) 


c -> 



VI. 




The final assumptions are that the total current density is a sum of a supercurrent 
density and a normal current density which is connected to the electric field 
through Ohm’s law. Thus, 

VII. 7 = 7s + Tn 

VIII. 7i ^ ■ 


2. The cgs Gaussian system of units will be used for calculations in this memorandum; key 
results will be expressed in the SI system of units. 

3. As in Lorenz’s electronic theory, lower case symbols are used to denote local values of 
internal fields. 
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These equations can be combined to give 


c . 9p 

— p + Atto — + 

at 


a^p 

'a?' 


0 


and 


ah 


a^h 


(curl curl h + — h) + 47ra — + — -- = 0 

x^ at at^ 


with equations for j and e identical to the one for h. For quasi-stationary conditions 
(that is, for frequencies oj satisfying the inequality cj«c^/4ttoX^ ^10*^ s"' ), 

div e = 47Tp - 0. The equation for h reduces to 


curl curl h + — h = 0 


with identical equations for j and e . 

To complete the theory, boundary conditions must be given for all surfaces of 
discontinuity where different bodies border each other and where the constants 
characteristic of the material change discontinuously. The boundary conditions are that 
h||, e||, and hj^,jj^ must be continuous; hj^ is continuous because the permeability p has 

been taken as unity. Finally, (X^jg)|j must be continuous at the boundary between two 
different superconductors. 

This brief review of the London theory was included to make this memorandum 
self-contained; a fuller presentation with sample calculations is given in the excellent 
monograph by F. London [3]. 


B. Rotating Superconducting Sphere 

For the case of a rotating superconductor, it is assumed that equations (V) and 

(VI) for the supercurrent are still valid. Using tlie relations X^ = m*c^ /dTro^e*^ and 
— ^ ^ 
jg = e*ngVg, these equations can be rewritten in terms of the velocity field Vg of the 

superelectrons: 
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curl Vg = -(e*/m*c) h 


— = (e*/m*) e 


where e* = 2e and m* = 2m are the charge and mass, respectively, of a superelectron 
(Cooper pair). The charge density p and the total current density j are 

p = e*(ng-ngo) 

T = e*(ngVg - ngQVQ) 


where n^ is the density of superelectrons -in the rotating superconductor, n^Q is the 

density in the same superconductor when stationary, and vq is the local state of motion 
of the body. 

Consider now the rotating superconducting sphere. To balance the centrifugal 
force on the electrons, the internal electric field must be 


->• mco -> 

e = (x + y) 


Note that p = (l/47r)div e = (mo;^/27re) = (m*co^/27re*), so that 


"s “ "sO 
"sO 


m*co2 



10 ’ 3 ' 


Thus (ng-ngQ)/ngQ is much less than 1 even for very large co, and to an excellent 
approximation 

T = e*ngo(Vg - Vq) . 
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Then 


curlj = e*ngQ(curl Vg - curl vq) = e*ngQ(-e*/m*c)h - e*ngQ(2co) 


Neglecting the displacement current in the Maxwell equation (I) and operating through 
with the curl, one has (using the above); 



Also, 

curl curl j = -e*ngQ(e*/m*c) curl h = -(47rngQe*^/m*c^) j 


so that' 




curl curl j 



Since boundary conditions are the same for the rotating superconducting sphere as for 
the stationary sphere, one can now prove a simple theorem which will be extremely 
useful for torque calculations: 

Theorem: 

— >■ — 

Let h^ and jr» be the solutions to London s equations for a stationary 

superconducting sphere in an arbitrary field B, and let hj^ and be solutions for the 

rotating sphere in zero applied field. Then the solutions to the London equations for the 

— ^ ^ ^ 

sphere rotating in the field B are h = hQ + h^ and j = Jq Jl- 
Proof : 


X^ curl curl j = X^ curl curl Jq + X^ curl curl = -Jq “ Jl “ “J 
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- ^ y y / ^ 2m ^ c ^ 

curl curl h = X^ curl curl 1 iq + X^ curl curl hj^ = - 1 iq -(11^ + — ojj 


-> 2m*c 
= - h : — io 


It is clear from symmetry that the theorem also holds in the case of a spherical 
superconducting shell. 

The current and field solutions for the rotating superconducting sphere, and for 
the stationary sphere in a uniform'* applied field, are derived in London’s book [3]. For 

the stationary sphere in a uniform applied magnetic induction Bq, one has (with the 
z-axis parallel to Bq as in Figure 1): 


h ° = 
“r 




- 2M' \ 

^0 + — 


' M' \ 

- ^0 TT ) 


h/ = 0 


r > R , 


2A' 

h ° = — — (sinh Pr - Pr cosh |3r) cos0 
P‘‘ r 


A' 

~ [(1 + |3^r^) sinh Pr - Pr cosh j3r] sin0 

j3 r^ 


r < R 


h/ = 0 


4. Because of the diamagnetism of the superconductor, h 9^= Bq outside the sphere. 
However, it is required that h ->■ Bq as |r| ->• <». 
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where 


|3^ = 1/X^ = 47t nge*^/m*c^ 

M, , . , 3 .cosh_gR ^ ^ \ 

2 y |3R sinh /3R R^ / 


2 sinh |3R 

— ^ 

For the rotating sphere, one has (with tire z-axis chosen along co as in Figure 2) 


L 2M 

hr = — COS0 


L M 

H = p-sin0 ) 


h '~ = » 


r > R 


L _ m’^c 
^r - — 


2A , 

2(0 + — ^ (sinh |3r - ]3r cosh j3r) 


COS0 


^0 


L m*c 


?[< 


2cj + — I ( 1 + r^ ) sinh j3r - (3r cosh jSr J | sin0 


r < R 


h^*^ = 0 
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MAGNETIC FIELD 
\ LINES 



Figure 2. Coordinate system for a rotating superconducting sphere witli no field applied. 
(The X axis of the coordinate system points out of the page in this figure.) 


and 


JL 


3 nge*coR j / I 

— I cosh Sr sinh Br 1 sin0 

sinh |3R |3r \ j3r 


A 

•(p 


where 


3cuR 


j3^ sinh |3R 


and 


M 


m*c 




3 cosh j3R ^ 3 

^R sinh (3R 


(M is called the London Moment.) In both cases the origin of the coordinate system is in 
the center of the sphere. These solutions will be used in ensuing calculations. 


C. Torque on a Rotating Superconducting Sphere in a Uniform 

Applied Field 

Although the rotor will be enclosed in a superconducting shield, the shield will 
have holes for electrode leads, pump out lines, gas spin-up ports, etc.; magnetic field 
could possibly enter through these holes. To obtain an estimate of the torque due to 
leakage field, the torque on a rotating superconducting sphere in a uniform applied field 

Bq will be calculated. This problem can be solved exactly. 


The torque N on a current distribution j(r) in a magnetic flux density h(r) is given 
by [10]: 


N = - J d^r [r X (j X h)] 


^ ^ ^ ^ ^ 

Appealing to the theorem, one writes j = Jq + Jl and h = 1iq + h^, where j'q, jj^, 1iq, 
and h^ are “the solutions quoted at the end of the last section. Substituting, one has 


N 


1 

c 


f dh [7x (J) X ho)] + ^ 

V 


f d^r [7x (jl X ^l)] 

V 


+ 


1 

c 


/ d^r [r’x ^ X hL)] 

V 


-t- 


1 

c 


f d^r [7x 0^ X ho)] 

V 



It is easy to show by straightforward substitution that the first two integrals are zero, so 
that 


N = - / d^r [r X Go .X hj^)] + ^ / d^r [r x Gl x hQ)] = 

V V 


In a coordinate system with the z-axis chosen parallel to Bq, Jq has the form (see 
Section IIIB) 


-> A 

JQ ^ ^ 


and ho has the form (r < R) 


hQ = hj.° (r, d)r + h^“ (r, 6)d = Fo(r) cosd r + Go(r) sin0 6 


Similarly, in a coordinate system with the z-axis chosen along cj, has the form 


j’l = ^ 


and hj^ has the form 

% 

h^ = hj.^ (r, 0)r + hg^ (r, 0) 0 == cos0 r + G^Cr) sin0 0 

Now torque is a pseudovector and, hence, is invariant under coordinate rotations. 
Therefore, in evaluating the integrals, the orientation of the coordinate system may be 

chosen for mathematical convenience. In the integral choose z parallel to cv with y 

perpendicular to the plane defined by co and Bq; in the integral choose z parallel to 
—►a 

Bq with y again perpendicular to the plane defined by cu and Bq. (After calculations are 

complete, the results must be expressed in the same coordinate system before being 
added.) 
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Thus in each integral the field has the form previously indicated, and the current 

— >■ — ^ 

must be rewritten in the new coordinates. If is the aiigle between cu and Bq, then in 

either case the transformation represents a rotation through angle CL about the y-axis: 
Clockwise (+J2) for one integral and counterclockwise (-J2) for the other. In Cartesian 

coordinates, one has for either integral (since r = in is invariant under rotation) 
h = sin0 cosd coscj) [F(r) + G(r)] x + sin0 cos0 sin0 [F(r) + G(r)] y 
+ [cos^0 F(r) - sin^d G(r)] z = x h^^ + y hy + z h^ 


and 


— > A A 

j = - f(r) sind sin0 cosH x + f(r) [sin6 coscj) cosCl - cosd sinf2] y 


A A A A 

+ f(r) sin0 sin<^ sinf2 z = x + y jy + z 


(Since the two integrals have the same mathematical form, the subscripts will be 

— ^ ^ ^ ^ 

temporarily suppressed, with the pairings '••nderstood.) Then for 

either integral, 


n( 1) = X - J d^r [r sin0 sin0 hy - r sin0 sin0 jy h^^ - r cosd Jz + r cosd h^] 


A 1 r 

+ y- j d^r [- r sin0 COS0 hy + r sin0 COS0 j h^ + r COS0 jy h^ - r COS0 h ] 


A 1 p 

+ z- Jd^r[rsin0 cos0 j^ h^ - r sin0 cos(^ j^^ h^ - r sin0 sin0 jy h 


y “z 


+ r sin0 sin0 j^ hy] 

I 
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— ^ ^ 

Substituting for the components of j and h from above and integrating over d and 

(f) gives 


= 


A 

y 


47T 

3c 


sinfi 


R 

f dr r^ f(r) F(r) 
.0 


Thus each of the torques lies along the y-axis of its coordinate system. Since ihe 
y-axes coincide, the two torques can be added (after the change for the second 

integral) to give the total torque: 

= ^( 1 ) + 1 ^( 2 ) = 5 


R 

f r^ fj^Cr) pQfr) dr 
0 


47T 

+ — sinJ2 
3c 


47T 


R 


— sini^ J" r^ fQ(r) F^fr) dr 


The term can be viewed as the shielding current interacting with the London field 

and the term as the London current interacting with the external field penetrating 
the superconductor. 

Substituting from Section IIIB yields 


nO) = 


^ to Bq sinL2 


R 


sinh ^R 


R 

J r(sinh j3r - jSr cosh jSr) dr 
0 


3R^ r dr 

,2 j p- ^r - Pt cosh /3r)^ 


|3^ sinh^ |3R 


0 


and 


n(2) 


m*c 


cj Bn sinLZ 


3R^ 


R 


dr 


— u; un 

^ U Q2 Cinh2 


|3^ sinlr 13R 


f — (sinh |Sr - j3r cosh j3r)^ 

J j-2 
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Thus, 


INI = 


m*c 

— — oj Bq sinfi 


R 

sinh |3R 


R- 

y r(sinh j3r - |3r cosh |3r) dr 

0 


Carrying out the integration, one obtains the exact result: 


■ INI 


nv^c 


R^ Cl) Bq sinS7 


/ 3 cosh |3R 

y j3R sinh j3R 



The result can be written in vector form as: 



For a rotation frequency of 200 Hz (co = 4007 t rad/s), a field strength of 1 
microgauss, and a rotor radius of 2 cm, the maximum torque is 


N 


max 


e'^ 


R^ Cl) Bq — 5 X 10'*° dyne-cm 


Nmax ^ 5 x 10-’ Nm 


Note that 


X 


5 X 10'^ 
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which means that the interaction of the London current witli the external field 
penetrating the superconductor is negligible in comparison with shielding currents 

interacting with the London field. This is quite reasonable: In and Iiq both 

become exponentially small at distances from the surface which are large compared to X, 
but in only Jq becomes exponentially small. 


D. Rotating Superconducting Spherical Shell 

The solution to London’s equations for a hollow superconducting sphere rotating 
with constant angular velocity (no external fields) has been calculated by J. H. 
Derrickson^ of NASA and others^ [11]. The results are: 


j 


s 

[A(sinh |3r - jSr cosh j3r) 
r 


+ 


B(cosh j3r - pr sinh j3r)] sin0 


A 

0 


2M 

hj. = -y COSd 


M 

hff = —sine 




m*c r 2A 
1 = 2co + ^ 

r e* L 


r > Rr 


2B 

(sinh |3r - j3r cosh j3r) + — (cosh |3r - |3r sinh |3r) cosO 


A 

? 


h^ = — — I - 2co + -J [(1 + |3^r^) sinh |3r - |3r cosh j3r] 


+ -r [(1 + |3^r^) cosh |3r - |3r sinh |3r] ) sin0 


Ri<r<Ro , 


h0 = 0 


5. Derrickson, J. H.: The Magnetic Field Generated by a Superconducting Spherical Shell 
Rotating with Constant Angular Velocity. MSFC IN-SSL-N-70-9, July 27, 1970. 

6. Pondrom, W. L., Jr.: London Moment of a Rotating Hollow Superconducting Sphere. 
Unpublished Report. 
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r < Rj 


hj. = Hq COS0 
li 0 = - Hq sin0 
= 0 

Here Rq is tlie outer radius of the shell, Rj is the inner radius and the constants M, A, B, 
and Hq are given by: 


3cuRq ■ [(3 + ( 5 ^ Rj2 ) cosh jSRj - 3 jSRj sinh |3Rj] 

~~ F ~ [(3 + |32 Rj^) sinh /3(Rq - Rj) + 3 )3Rj cosh j3(RQ - Rj)] 

3gjRq [3 j3Rj cosh ^Rj - (3 + i3^Rj^) sinh |3Rj] 

|3" [(3 + |3^ Rj") sinh i3(Rq - Rj) + 3 jSRj cosh i3(Rq - Rj)] 


2 m*cco ^ ^ ^^0 

^ ~ (3 + sinh |3(Rq - Rj) + 3 /3Rj cosh /3(Rq - Rj) 


m*c wRq^ 

3 3 

1 + 

(3+^^R7)cosh|3(RQ-Rj) + 3^Rj 

sinh /3 (Rq - Rj) 

1 

e* 

P'Rq" P^o 

(3 + ]8^ R3 ) sinh |3 (Rq - Rj) + 3 (]R. 

cosh.)3(RQ - Rj) 

j 


In the limit Rj 0, these results approach those computed by London for a rotating 
solid sphere of radius Rq (see Section IIIB). In the limit j3Rj » 1 and |S(Rq - Rj) » 1 

the result approaches that for a solid sphere whose radius is much larger than the 
penetration depth. 

A limit of special interest is /3Rq» 1 with Rj ->• Rq. Calling Rq - Rj = d, one has 


Hq 


2 m*ccj 


Rq d/3) 

1 + (P^ Rq d/3) 


2 m*ccj ^ ^0 
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m*c , 
M — cjRn^ 


Rq d/3) 

1 + (p^ Rq d/3) 


m*ca)RQ^ P^ Rq d 


e"^ 


ns 


coRq sin0 <p 


In this limit the current density and fields approach those of a rotating charge 
shell of density a = -nse*d and total ch- rge Q = 47rRQ^dnse*. The superelectrons are thus 

nearly stationary in the lab frame and move with a velocity Vg = -cuRq sin 6 with respect 

to the lattice, so that a reasonable requirement is that the current density not exceed the 
critical current density; i.e., it is required that |nse*coRQ| < where is the critical 

current density. From the Landau-Ginzburg theory for a thin plane superconductor, one 
obtains [4] ; 


Jc(T) 


C $Q 

12v^7T^ \HT) ^(T) 


where 4 >q is the flux quantum. Using nse*coRQ = m*c^coRQ/47re*X^ , the requirement 
becomes 


f 


CO 

— < 
2tt 


mcRQ ^(T) 


0.018 

^(T) 


Hz 


Therefore, it appears that a thin rotating shell gives the full London moment, provided 
the rotation frequency is not too large and provided the pair density is essentially that of 
a bulk superconductor. The latter condition is equivalent to the requirement that 


d > 1(0) ~ 400 A 


E. Effect of Excess Charge on the Rotor 

An excess charge on a superconducting sphere at rest will be uniformly 
distributed on the surface. The requirement that the tangential component of "e be 
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continuous at the surface means this will also be true for a rotating superconducting 
sphere. If the total charge is Q, the surface charge density is a = Q/47 tR^ . At steady state 
this charge rotates with the lattice, constituting a current density 

A 

J = a r 60 sin0 5(r - R) 0 


The vector potential A is then 


Xf, = i f,S , 

ir - r 


The expansion of ir - r f in spherical harmonics is { 10] : 


r - r 


oo g r ^ 


e=o m=-e 


2C+1 jC+l 


— >■ — A A 

where r^ (r^) is the smaller (larger) of |r| and Ir |. Using the relation 0 = - sin0 x + cos0 




the evaluation of the integral for A(r) is straightforward. The result is 


47t 060 R'* a 

A(r) = sind 0 , i > R 

c 3r 


and 


A(r) 


47t 060 Rr A 

sind 0 

c 3 


r < R 


From B = curl A the 


magnetic field is 
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< 



B, 


B. 


2Qcj 

3rT 

Qco 



COS0 


sin0 


- 0 


B, 


2Qco 

3Rc 


COS0 


B0 


2Qcj 

3Rc 


sin0 


B. 


0 


r > R 


r < R 


The field inside the charge shell is a uniform field parallel to co; since this field is 
due to “external” sources, the solutions to the London equations are those of a sphere 
rotating in a uniform applied field. By appeal to the theorem, one can then write down 
the solution from the results of Section IlIB: 

For r < R, 


r g* 


m^c 

^0 “ 


2co + 


6cjR 


(Qe*/R) 


I 


|3^r^ sinh |3R 3m*c^ J 


1 - - — r— : — I (sinh j3r - jSr cosh j3r) 


COS0 


- 2cu + 


3cjR 


sinh jSR 


I 

! 


(Qe*/R) , 

1 - — — 1 + j3^ r^ ) sinh |3r - j3r cosh j3r] 


sin0 


h0 = 0 


The field outside the rotating superconductor will be the superposition of the 
field due to the supercurrent and the field due to the rotating charge layer: 
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h 


r 


2m*ccj 


R \^ / 3 cosh PR 3 \ 

7/ V ~|3Rsinh/3R ^ p^J 


cosd 


3Rc \ r/ \j3Rsinh/3R p^Ry 


h0 


m*ccj / R V 
e* \7 / 


3 cosh j3R 
i3R sinh |3R 



sin0 


gQ /rV /3coshgR , 

3Rc y.r y \^j3Rsinh/3R jS^R^I 


h0 = 0 


For Q = Ne the relative magnitudes are 



m*c , 

cuR 3 
e* 


^ mc^ R 1 

. io‘" 

Qcu R^ 3 cosh j3R 
3Rc (3Rsinh|3R 

(e7R) X N 

N 


so the direct effect of a small excess charge is negligible. These results are easy to 
understand; The rotating surface charge generates a dipole field outside and a uniform 
field inside the superconductor; a supercurrent flows to exclude (cancel) the uniform 
magnetic field from the interior of the superconductor, generating a uniform field inside 
and a dipole field outside. The net result is that these fields cancel everywhere to at least 
order X/R. 

For completeness, the current density j in the superconductor is 


3 n^e* cjR 
P^Y^ sinh(3R 


Qe*/R\ 

3m*c^y 


A 

(sinh Pi - jSr cosh |3r) sin0 ' <p 


Josephson has argued [12] that m* is the rest mass corrected by the work 
function W, the energy difference between free space and the Fermi surface; i.e., that it 
is the total energy, m*c^ = 2mQC^ - 2W, necessary to create a pair in the metal which 

determines the relevant mass m*. Since the addition of charge Q to an isolated 
superconducting sphere increases the potential by Q/R, the work function becomes 
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W = Wq + Qe*/2R. If Josephson’s argument is correct, this means the addition of charge 
could significantly change the London moment through m*. 

IV. CALCULATIONS USING CLASSICAL ELECTRODYNAMICS WITH THE 
ASSUMPTION OF A PERFECT MEISSNER EFFECT 


A. Introduction to a Macroscopic Approach 

The Meissner-Ochsenfeld experiment indicated that the magnetic induction B 
vanishes inside an ideal superconductor. Meissner’s result is contained in the London 
theory, with the plausible restriction that the magnetic induction does not vanish 
discontinuously at the surface, but instead it decreases to nearly zero over a distance of 
order A, the penetration depth. The tjiemodynamic magnetic induction ^ is defined as 
the volume average of the local field so that for large bodies, ^ = 0 in the interior. 

— >■ 

The general boundary condition for the magnetic induction B at a surface is 


— > — > A 

(B 2 “ B j ) ■ n — 0 


A 

where n is a unit vector normal to the surface. If the Meissner result is used, then, with B 
the external induction and n the outward normal, the boundary condition at the surface 

A 

of a superconductor becomes B * n = 0. Use of this boundary condition for 
superconductors whose dimensions R are large compared to X gives values for the 
external magnetic induction which are correct to order X/R. 

— ^ 

As an example, consider the case of the sphere in uniform applied field Bq. The 

external magnetic induction can be derived from a scalar potential T>. Choosing the origin 

— 

of the coordinate system at the center of the sphere with the z-axis parallel to Bq, 
symmetry requires that <I> have the form 


T> 


Yj r*^ + j3g r Pg(cos0) 

e-0 ^ 


where the functions Pg(cos0) are the Legendre polynomials [10]. As |r|^°°, 


B - - grad <h 


A 

- Z - 


grad 


Y Pg(cos0) 

e=2 
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but as ir|->oo, B Bq, which requires that - -Bq and ag = 0 for Ol.Thus, $ 
becomes 


fh = - Bq r COS0 + ^ r Pg(cos0) 

C=0 


At the surface of the sphere, B • n = Bj.(r = R) = 0, so that 


0 


3<l> 



Bq COS0 + X' (C+1) i3g p^(cos0) 

C=0 


Orthogonality of the Legendre polynomials gives iSg = 0 for £ ^ 1, and (3j = -(BqR^/ 2). 
Hence the result 



It is seen that this result agrees with the exact result (Section lllB) to order A/R. 

— > A 

The condition B ■ n = 0 holds only for externally applied fields in the case of a 
rotating superconductor; the rotating superconductor generates a magnetic induction 
which has a nonzero normal component. For a rotating superconducting sphere with the 
origin of the coordinate system at the center of the sphere, the boundary condition for S' 
becomes 


-> A 2M A A 
B • n = — (co • r) 
R"* 
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where M = (m* c/e*)coR^ . This reduces to B ■ n = 0 for co = 0. Inside the rotating sphere, 
the tlrermodynamic magnetic induction is 


2 m*c ^ 
B = : — CO 


which reduces to B = 0 for co = 0. 

In the complete Meissner region, tlie thermodynamic magnetic field H is 

undefined. One may define H [and concomitantly the magnetization m through 
— ^ ^ 

47 rm = (B-H)] so long as the definition , is consistent with thermodynamics and 
macroscopic electromagnetic theory. From Maxwell’s equation (I) at steady state, one has 
for any surface S bounded by a curve C, 


r -> -i- 47t f ^ ^ 

9 H • dC = — J j • ds 
C ^ S 


in a superconductor, currents flow within order X of the surface, so that j = 0 deep within 
the superconductor. Then for any arbitrary small closed path C deep within the 
/* — ^ ^ ^ ^ 

superconductor, 9 H • dfi = 0, and the choice H = B = 0 in the interior is self consistent. 
C 

The general boundary condition for H is 


A 

n X 


^ 47T ->■ 

(H 2 -H 1 ) = -K 


where K is the surface current density. The boundary condition thus becomes 


A 

n X H 


47T -> 
— K 
c 


at the surface of a superconductor, where H is the external magnetic field. 


24 



Continuing with the example of the sphere, one has H = B (p = 1) outside the 
superconductor, so that 




c A 

— Bg (r = R) 0 
47T 


or 



From Section IIIB, one has for the current density 


3B 


0 


R 


J - 


47t 2r^ sinh |3R 


(sinh ^r - )3r cosh |8r) sin0 0 


Integrating 


R 


/ J0dr 
0 


c R sin0 

47t 2 sinh ^3R 


R 

J (r‘^ sinh j3r - |3r'* cosh j3r) dr 
0 


c 

4tt 



sin0 


j3R \ , c 3^0 

sinh/3R / » 1 iff T~ 


The current density j is appreciable only within a distance of order X of the surface, and 

— >■ C A 

tlie equation K = — [nxHl implies that this current is treated as a true surface current. 
4ff ^ 

That is, to within X/R, K is just the integrated value of the current density obtained by 
solving London’s equations. (The apparent difference in sign is caused by tlie fact that j 
is the electron current density 0 e* < 0), whereas the macroscopic theory, as usual, 
assumes positive charge carriers.) 
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For the rotating sphere it will also be assumed tliat p = 1 (i.e., take B 
interior, so that the boundary condition becomes 


H) in the 


n X [H - (2m*c/e*) co] 


(47t/c) K 


where H is the external field. The choice p = 1 is consistent with the London theory as 
well as thermodynamics and classical electrodynamics. It is seen that the stationary 
superconductor (co = 0) is a special case of the above. 

The mode of attack for torque calculations will be as follows; The magnetic 
induction B in the region between the rotor and the shield will be calculated as the 
gradient of a scalar potential and will be evaluated using the boundary conditions; 


-> A 2M A A 
B • n = — (cu • r) 


at the surface of the rotor, 


B • n = 0 at other superconducting surfaces. 


Since B = H, the rotor surface currents can be calculated when B is known; 


K = — 
4ir 


A 

n 


/ 2m*c 


j 

The torque can then be calculated using 


N 


R3 

c 


2tt 7T 

f d<P f ^ ^ 

<^>=0 0=0 


As a final example, the torque on a superconducting sphere in a uniform applied 
A — ^ A 

field is calculated. Choose z parallel to Bq and y perpendicular to the plane defined by 
— ^ ^ ^ 

Bq and co. Requiring B ^ Bq as r ^ and using the boundary condition 
B • n = (2M/R^) (w '^) at the surface of the sphere, the scalar potential is 


26 



BoR= 


47t M 


<I> = ~ Bq COS0 + — — ^ 


2r^ 


T7I E Yim(0,0) 


m=-l 


BqR^ 


M 


- - Bq r COS0 - - - - 1- cosd + — (sinS2 sin0 cos0 + cosfi cos0) 
2r 


where is the angle between Bq and co, and tlie functions Y^^id,<j>) are the spherical 
harmonics as defined by Jackson [10]. 

The components of the magnetic induction B = -grad are 


■■■('-a 


r. . 2M 

Bq COSO + (sinJ2 sin 6 cos 0 •+ cosfZ cos0) 


B, 




T. • . M 

Bq smd - — (sini2 cos0 cos0 - cosf2 sin0) 


M 

B 0 = — sinf 2 sin 0 


In vector notation, 


B 


Bnz - 


BqR3 

2r3 


[3(P' 


A A A 

z) r - z 


] + 


M A 
r-* 


A A A 

r) r - cu] 


The internal magnetic field H is H = (2 m*c/e*)cu 


2M A 

- CO, so that the surface current is 


K 


C A 
— r X 

4ir 


Be 


^0 A 

-- I3(r 


A A A M A A;_ A A 2M A 
z) r - z] + -^ [3(r • cu) r - co] - — r w 


c 

47T 


A A 3M A A 

— - (r X z) - — (r X co) 

2 R"’ 
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The torque N is then given by 


N = 


R3 

c 


2it 

/ 

0=0 


7T 

^ A 


'3Bo 

A A 

3M A A 

d0 J dd sin0 r x ' 
0=0 

47T 

2 

(r X z) - 



A^O AAAA M AAAA 

Bn z - — [3(r • z) r - z] + :^ [3(r • cj) r - cu] 
'^2 R 


There are six vector products in N: 


A AAA AAAA 

1. r X [(r X z) X z] = (r ‘ z) (r x z) 


A AAA A 

2. r X [(r X z) X r] = (r x z) 


A AAA AAAA 

3. r X [(r X z) X co] = (cu • r) (r x z) 


A AA a AA AA 

4. r X [(r X co) X z] = (r • z) (r x co) 


A AAA A A 

5. r X [(r X co) X r] = (r x co) 


A AAA 

6. r X [(r X co) X co] 


AAAA 

(co • r) (r X co) 


Terms 1, 2, 5, and 6 integrate to zero. Integrating 3 gives 


2n 7T 

f - AAAA 

J -dip j dd sin0 (co • r) (r x z) 
0 0 


2rr n 

A z' z' AAA 

• z X dp I dd sin0 (sinfi sin0 cosp + cosS2 cos0) (x sin© cos0 + y sin0 sin0 + z cos0) 

0 0 


A 

- Z X 


2tt 


r t r A A 

sin^2 j d0 sin^0 j d0 cos0 (x cos0 + y sin0) 

0 0 


4n A 

= sin^2 y 

3 
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Similarly, integrating 4 gives 


2n 

f 

0 


Thus the torque 


/ 


dd 


A 

(z 


A A 
r) (r X 


A 

CO) = 


47T a 
— sinJ2 y 
3 


0 

N is 




R3 c 

c 47 t 2 




/ 3M\ 

/3Bq> 

\ 47T 

+ t) 


— sinS2 

\ rV 

V 2 ) 

3 


= - Bq M sinfl y 


In vector notation, 


N = M X Br 


M = (m*c/e*) cj 


which agrees with the exact solution to order X/R. 


B. Torque on a Rotating Superconducting Sphere in 
an Arbitrary Magnetic Field 

Using this macroscopic theory, the rotor torque in an arbitrary magnetic field can 
be calculated. Let S be the distance between tlie center of the rotor and the nearest 
external current or superconducting surface. Then in the volume defined by R < lr| < S, 
the magnetic induction B can be derived from a scalar potential 4> through B == - grad 4>. 
The boundary condition 


A 

B • n 


2M A A 


at the surface of the rotor requires the scalar potential to have the general form 
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<|.(r, 0 q) ^ Z YjJ (n,0Q) Yj,,(0,0) 

k=-l ■ 


+ 


°° £ r 

Z Z "£m 1 

(?i) 

£=1 m=-£ ■- 



Ygni(6,0) 


where is the angle between a; and the z-axis, 0 q is the angle between the x-axis and 

the projection of co in the x-y plane, M is the London moment [M = (m* c/e*)coR^], 
and the functions Y^^{d,(p) are the spherical harmonics. 

The scalar potential for an isolated rotating sphere is 


1 

Hl = yMr-^ E Ylk Yik(«.« 

k=-l 

— y — > 

so that one can write <J> = + <|)q. The magnetic induction B is B = -grad<l> = -grad^^ 

— ^ ^ ^ ^ ^ 

-grad<J>Q = Bj^ + Bq. (Bq = B - Bj^ should not be confused with the uniform applied field 
considered earlier.) The components of B are 

1 

94> Stt M • 

Br = -5- = T? E V,J(n,«.o)Y, ,<(«,«.) 

k=-l 


00 £ 

£=1 m=-£ 


10) , 


1 


1 3<I) 47t M V * 9 

k=-l 


00 £ 


Z Z "£m 

£=1 m=-6 


rC-1 + f JL )r 2£+1 r-(fi+2) 
£+1 
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1 ^ 

r sin0 00 


47tM 

3 sinfl 


k=-l 


oo 2 


- z z 


sind 



0 

00 


Ygi„(0,0) 


Since Bq is real, = (-1)^ o:g^_i 


nr 


Because B = Bq + Bj^, K = Kq + K^, where 


Kr 


C A -»■ 

— n X Bf 
47T ^ 


Kt = — 


C A 
n X 
47T 


B, 


2 m* c 


cu 


C 3M A A 
47t 


Quite generally, 


->■ a-^->^ k"a.|a ->■ ->• -> ->i 

N = — f dn' [r X (K X B)] = — j dS2 ^ r x [(Kq x B^) + (K^ x BQ)]j 


R 

c 


A — ^ ^ C 3 ^ AAA 

where d^2' is the differential of solid angle. But r x (K^ x Bq) “ ^ ’ r) (r >< <^) ~ 0> 

— ^ A • 

since Bq ' r = 0 at r = R. Thus 


N 



(Kq X Bl)1 


Since 


A 

r X 


(Kq X Bl) 


47T 


[d (- B, 


0 „L 


Bp + 0 {Bq 


0 


Bp] 
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N 


47T 


2iT 7T 

f d<P f dd Sind [0 (- bJ') + i (B^ bJ')] 
0 0 


The evaluation of the integral is straightforward but tedious. Recall that 


A A 
r • X 

= sin0 COS0 

A A 

0 • X = COS0 COS0 

A A 

0 ■ X = - sin0 

A A 

r ■ y 

= sin0 sin0 

A A 

0 ■ y = COS0 sin0 

A A 

0 • y = COS0 

A A 

r • z 

= COS0 

A A 

6 ■ z - - sin0 

A A 

0 • z = 0 



-> -> 



Substituting for the components of B^ and Bq, one has in rectangular coordinates: 


N., = 


2M 


oo g 


Z Z Z ‘ “em ^ik (^.0 q) / 

<=-l £=1 m=-e ' ^ 


= M i /— sin n (aj j - a|j) 


for 0Q = 0 ; 


N. 


k e m ' \ 

+ /d^'sin0 Yjk 


, COS0 COS0 
sin0 


Ylk^en, 



Nx = M i /— cos n (aji -ajj) 


for 0Q = 0 
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k fi m ^ ^ ^ 


COS0 cos<j> 
sin0 


Jdn'cos0 Yjj^ ^Yflrn^ 


Ny - M cos (aj j + aj j) + M Ojq sin 12 for 0 q = 0 


Using this general result the previous calculation can be verified. For a uniform applied 
field, 


4>0 




Yio(0>0) 


so that 



^Cm ~ ^ S. ^ I, m ¥= 0 . Then 


Nx ^ 0 , 

Ny = - M Bq sin 12 
Nz = 0 > 


as before. 

One additional point should be mentioned. Because a superconductor is 
diamagnetic, it experiences a force when immersed in a nonuniform magnetic field. This 
force can cause a torque if the center of support is not the center of mass. Using the 

— > A 

boundary condition B • n = 0 at the surface, tlie force on a superconducting sphere in an 
arbitrary magnetic field has been calculated by Harding [13]. 
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C. The Infinite Superconducting Plane 


Assume that the half-space defined by the equation x < 0 is occupied by a 
semi-infinite superconducting slab and that tliere are no other superconductors in the 

half-space x > 0. One asks for the induction B due to the infinite superconducting plane 
and an arbitrary current distribution in the region x > 0. 

A 

Since tlie boundary condition is B • n = = 0 at x = 0 and since there are no 

other superconductors present, this problem is a simple image problem. The solution is 
found by replacing the superconductor by the image of the current distribution mirrored 

in the y-z plane. The solution for B in the region x > 0 is then the superposition of the 
fields due to the original current distribution and the image current distribution. 


D. Torque due to a Current in the Readout Loop 

Choose a coordinate system with the origin at the center of the sphere and with 
z-axis perpendicular to the plane containing the readout loop. The center of the sphere is 
assumed to lie in this plane. Let the radius of the readout loop be R -i-e, where R is the 
radius of the superconducting sphere. The loop itself is assumed to have negligible cross 
section so that the current density of the loop is 


1 5(cos0) 5(r - R - e) ^ 

(b 

R 


where I is the current in the readout loop and 5(x) is the Dirac delta function. 

Assume for the moment that the rotor is stationary. Shielding currents flow in 
the rotor“'to exclude the field of the current loop. In the limit of tight coupling (e 0), 
the rotor appears locally to be an infinite superconducting plane. In other words, to have 
a perfect Meissner effect, the field of the supercurrent must cancel the field of the 

— y 

readout loop, so that in the limit e ^ 0, the supercurrent density Jq is just 


^ I 5(cos0) 5(r - R)a I 5(cos0) 5(r - R) a a 

Jq r 0 = [- sm0 X -I- COS0 y] 


Appealing to the theorem, for the rotating sphere, one has 


-> 1 5(cos0) 5(r - R) a a 

j “ Jl ~ [“ sin0 X -I- COS0 y] 

h - hL . 
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In the coordinate system with z' parallel to co (see Figure 3), the London field has 
the form 


hLC?) = r'F(r') cos0' + ^ G(r') sin0' 


- x’[F + G] cosd^ sind’ cos0* + y’[F-+ G] cos0’ sin0’ sin0^ 


+ z' j [F + G] cos^e' - G j 


Denoting the angle between z and z' by and choosing y and y' to coincide, the 
components of hj^fr ) can be rewritten in terms of the unpriraed coordinates: 



Figure 3. A current loop around the rotor. (The center of the sphere coincides with 
the center of the loop; the Y,Y' axes coincide and lie in the plane of the loop. 
The Z axis is perpendicular to the plane of the loop; the Z' axis 
lies along the spin axis of the rotor.) 
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(h^) = [F + G] cosH (sinJ2 sin0 cos0 + cosFi cos0)(cosJ2 sin0 cos<p - sinfi cos0) 


+ .[F + G] sinS7 (sinJ2 sin0 cos0 + cosJ2 cos0)^ - G sinfi 


(h]^) = [F + G] sin0 sin0 (sinJ2 sin0 cos<^> + cosfi cos0) , 


(hj^) = -[F + G] sinJ2 (sinfi sin0 cos0 + cosQ, cos0)(cosF2 sin0 cos0 - sinH cos0) 


+ [F + G] cos^2 (sin^2 sin0 cos0 + cosJ2 cos0)^ - G cosJ2 . 


The torque on the rotor is 


N = - j.d^ r[r X (j X h)j = - Jd^ r[r x Gq ^ h^)] 




Substituting the expressions for Jq and h|_^ into the above, carrying out the vector 
multiplication, and then integrating, one obtains the result 


N = - y - 7 tR 2 F(R) Sinn 
c 


Now 


2m*ccu , 

F(R) = + terms of order X/R , 

e* 


so that to lowest order in X/R the torque is 
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Thus to lowest order in X/R, the torque is the same as that on a current loop of radius R 
in a uniform magnetic field equal to the London field. Witlr I = 0.5 pA, R = 2 cm, and 
cj = 400 7T s'* , the maximum torque is 


Nmax ~ 9 X 10'*« Nm 


E. Torque due to a Flat Region in the Superconducting Shield 

A change in the gyro housing design was proposed in order to reduce difficulties 
in the precision machining of certain parts. However, the change in design necessitated 
departure from spherical symmetry for the superconducting magnetic shield, which would 
introduce extraneous torque on the rotor. In the new design, the superconducting sliield 
would be spherical except for a flat spot of half angle 0 q, as shown in Figure 4. 



Figure 4. A proposed shield design. (The shield is spherical except for a flat spot which 
subtends a half-angle 6 q at the center of the spherical portion. The 

spherical part of the shield is concentric with the rotor.) 
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The torque for this shield design with the spin axis at a small angle 12 with the 
respect to the axis of symmetry has been calculated by C. Ebner and C. C. Sung’. Their 
result is expressed in the form 


N = C(0q,R'/R) (m*cco/e*)’ R^ sinl2 


Here R is the ball radius and R' is the radius of the spherical section of the shield. 
(Actually the coefficient C depends on the angle 12 also, but is rather insensitive for small 
12.) Ebner and Sung provided numerical values of C(0 q,R 7R) for 15°<0g<45° and 
1.5<R'/R<3.0. 

To get a better physical feeling for the problem and to gain insight into the above 
result, an approximate calculation for the limiting case of a large flat spot was carried 
out. 


Since the London field falls off as a dipole outside the rotor, the currents in the 
flat part of the shield should fall off rapidly away from the center of the flat. Therefore 
for a large flat spot in a large shield (R « R' cos0q « R'), the torque should be 
approximately that on a point dipole near an infinite superconducting plane (Fig. 5), a 
problem easily solved by image methods. 



Figure 5. Image method for calculating torque due to the flat spot in the point 
dipole/infinite superconducting plane approximation. 


7. Karr, G. R. and Hendricks, J. B.: Fourth Quarterly Progress Report, Contract NAS8- 
29316, Appendix 1 — Magnetic Shielding Calculations, by C. C. Sung and C. Ebner. 
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The torque is given by 


N = M X B (image dipole) 


- M X 


3 n (^ ■ n) - M' 
|r - r '1^ 


8d^ 


= [3 cosJ2 sin(7r - fl) - sin(jr - 2^2)] y 


A 

sin y 


16d 


where y is a unit vector out of the page in Figure 5. Substituting d = R' cos0q and 
M = (m*coj/e*)R^ yields the approximation 



For R = 1.9 cm, R' = 4.3 cm, 6q = 57 deg, and oj = 400 tt Hz, one has 


N “ 1.2 X 10'*® sin2S2 Nm 


in good agreement with the Ebner-Sung result of N s 1.6 x 10'*^ sinll Nm for these 
values. 


Since tlie actual rotor/shield arrangement has the same symmetry as the point 
dipole/infinite plane, it can be concluded that the torque will be zero for = tt/2 as well 
as for S2 = 0 and S2 = tt, a result not foreseen from tlie small angle result of Ebner and 
Sung. This observation is important because, after the Ebner-Sung calculation had been 
completed, the design was changed again; in the new configuration it is planned to have 
the spin axis such that Q, — -njl rather than — 0. 
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F. Discussion 


As noted in the introduction, one of tlie design goals of the gyroscope relativity 
experiment is that the precession due to nonrelativistic torques be less than 1 ms of arc 
per year. This requires that the magnitude of the residual torque be less than 2x 10'*® 
Nm. 


The preceding calculations show that to meet this objective, the superconducting 
shielding must be designed so as to reduce the magnetic field in the interior to 
approximately lO'”* gauss (10"‘ * tesla) or less, and the field due to trapped flux must be 
equally small. These restrictions were noted earlier by C.W.F. Everitt [14,15]. It should 
be possible to eliminate trapped flux by careful cooling in a sufficiently low field; with 
special techniques, field levels as low as 6 x 10'® gauss have been obtained [ 16] . 

Since the magnetic flux through a hole in a superconductor cannot change, any 
magnetic field leaking through a hole in the superconducting shield surrounding the rotor 
will be strongly attenuated. A ballpark estimate can be obtained for leakage of magnetic 
field through a hole in a spherical shield from data for superconducting cylinders. For a 
cylinder of radius a, the axial field decreases as B = Bq exp(- 3.4 z/a) and the transverse 

field decreases as B = Bq exp(- 1.8 z/a), where z is measured from the end of the cylinder 

[17]. If one lets z = R'-R, the distance between rotor and shield, then for a transverse 
attenuation of 10'^ an estimate of a ~ (R'-R)/9 for the radius of the largest permissible 
hole is obtained. For R' = 4.3 cm and R= 1.9 cm, tliis gives 0.3 cm. Since each 

hole attenuates more or less independently and since tire torque calculation assumed a 

uniform field, the number of holes should not be critical. The exponential z-dependence 

above results from separation of variables in cylindrical coordinates and the vanishing of 

— >• 

the normal component of B at the inner cylinder wall; a different functional dependence 
for the variable r would certainly be expected for a hole in a spherical superconducting 
shell. Moreover, the attenuation should depend on the size of the hole relative to both 
the shield thickness and the shield radius, probably through multiplicative factors. 
Finally, holes remove spherical symmetry and, like the flat spot, produce a diamagnetic 
torque. Clearly a more precise calculation of the hole problem should be undertaken. 

The torque due to the flat spot will depend on the pointing accuracy. For 

= 0 + (or .= ^ + 5f2), one has for 0 q = 57 deg: N 2.8 x 10'*® N m for 5f2 = 1 

deg, N =“ 4.7 X 10'^\^N m for Sfl = 1 min, and N ^ 7.8 x 10'^^ N m for = 1 s. Thus 
for a pointing error of less than 1 min of arc, the torque due to the flat spot will be 
completely negligible. 


8. Hendricks, J. B.: A SQUID Readout System for a Superconducting Gyroscope. 1974 
Applied Superconductivity Conference, Sept. 30-Oct. 2, 1974, Oakbrook, 111. (Proc. to be 
published). 
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The magnitude of the torque due to a current in the readout loop is mildly 
disturbing. In the direct SQUID configuration® , a dc current as large as that 
corresponding to half a flux unit (<Tq = 2.07 x 10'*® weber) can exist in the loop. For a 

SQUID inductance of 1.5 x 10'^ henry, this current could be as large as 6.9 x 10'^ A, 
with a corresponding torque of 1.2 x 10"'’ N m. However, a slow roll of the spacecraft 
about the optical axis of the telescope is planned. If the spin axis of the gyro coincides 
with the optical axis of the telescope, the average torque due to a current in the read-out 
loop will be zero. 
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